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ABSTRACT 

We demonstrate a new statistical method of determining the global photometric properties of the 
Milky Way (MW) to an unprecedented degree of accuracy, allowing our Galaxy to be compared 
directly to objects measured in extragalactic surveys. Capitalizing on the high-quality imaging and 
spectroscopy dataset from the Sloan Digital Sky Survey (SDSS), we exploit the inherent dependence 
of galaxies’ luminosities and colors on their total stellar mass, M*, and star formation rate (SFR), 

M*, by selecting a sample of Milky Way analog galaxies designed to reproduce the best Galactic M* 
and M* measurements, including all measurement uncertainties. Making the Copernican assumption 
that the MW is not extraordinary amongst galaxies of similar stellar mass and SFR, we then analyze 
the photometric properties of this matched sample, constraining the characteristics of our Galaxy 
without suffering interference from interstellar dust. We explore a variety of potential systematic 
errors that could affect this method, and find that they are subdominant to random uncertainties. 

We present both SDSS ugriz absolute magnitudes and colors in both rest-frame z=0 and z=0.1 
passbands for the MW, which are in agreement with previous estimates but can have up to ^3x lower 
errors. We find the MW to have absolute magnitude — 5log h = —21.00^03® and integrated color 
0(g-r) = 0.682™, indicating that it may belong to the green-valley region in color-magnitude 
space and ranking it amongst the brightest and reddest of spiral galaxies. We also present new 
estimates of global stellar mass-to-light ratios for our Galaxy. This work will help relate our in-depth 
understanding of the Galaxy to studies of more distant objects. 

Subject headings: Galaxy: evolution — Galaxy: fundamental parameters — methods: statistical — 
stars: formation — Galaxy: stellar content 


1. INTRODUCTION 

Galaxy evolution studies primarily rely on observa¬ 
tional comparisons between objects in the local universe 
(e.g., the Mi lky Way; MW) and those at higher red- 


shift, z (s.g. , Fab er et al. | |2067 Ilbert et '^|2010 Leau- 


thaud et al. 2UI2). h'or rnost galaxies of known z, rest- 


frame colors and absolute magnitudes are some of the 
easiest global properties to measure, regardless of their 
distance from us. Consequentially, color-magnitude dia¬ 
grams (CMDs) provide a fundamental tool for interpret¬ 
ing galaxy evolution, especially at large z where mor¬ 
phological information is difficult to obtain. Yet to this 
day, the MW’s position on such a diagram has remained 
poorly determined, despite being the galaxy we can study 
in the most detail. Due to our location embedded within 
the disk of the Galaxy, interstellar dust obsc ures stars 
and h ides most of the MW from our view (cf. Herschel 
17851. Furthermore, because bluer light is absorbed and 


scattered more efficiently by the d ust, the optical col- 


1989 


ors of distant stars are altered (e.g., Cardelli et al. 
Schlegel et al.||1998D. Thus, determining the global opti¬ 
cal properties of the MW from direct photometric obser¬ 
vations has proven extremely difficult, requiring uncer¬ 
tain corrections and assumptions that are vulnerable to 
error. 

For this reason, the history of measurements of the 
Galaxy’s global photometric properties is sparse. |van| 
der Kruit (1986, hereafter vdK86) made the most re¬ 


cent significant measurement, utilizing a novel technique. 
The Zodiacal cloud, a thick disk-shaped concentration 
of dust lying in the ecliptic plane (or zodiac), produces 
a glow of diffuse optical light throughout the nig ht sky 
via the reflection of sunlight (Reach et al. |1996 1. This 
diffuse glow, known as the Zodiacal light, introduces a 
significant amount of contamination to attempts to es¬ 
timate the amount of starlight from the Galaxy. The 
Pioneer 10 spacecraft, launched in 1972 on a mission to 
Jupiter, became the first space probe to travel beyond 
the asteroid belt, to distances where the effect of the 
Zodiacal light becomes negligible, van der Kruit used 
photometric measurements of the Galactic background 
light in broad optical blue (3950-4850A) and red (5900- 
6900A) bands taken by Pioneer 10, corrected for diffuse 
Galactic light and extinction, and compared to stellar 
distribution models in order to find Mb = —20.3 ± 0.2 
and B — V = 0.83 ± 0.15 in the Johnson magnitude sys¬ 
tem. 

Two earlier studies used a model of the Galaxy that 
consisted of a disk and spheroid component, but utiliz¬ 
ing different data and assumptions, in order to infer the 
luminosity and color of the Galaxy; both of these yielded 
significantly bluer color estimat es for the MW. First, 


de Vaucouleurs & Pence (1978 hereafter dV&P) had 


used a two-component model constrained to match the 
observed distribution of globular clusters in the Galac¬ 
tic bulge and the star counts near the Galactic poles 
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of the disk in the solar neighborhood in order to infer 
B—V = 0.53±0.05. This work also yielded Mb = —20.08 
(varying the shape of the bulge within this model yielded 
values ranging from -20.04 to -20.12); however, this es¬ 
timate was upd ated to Mb = —20.2 ± 0.15 in |de Vau^ 


couleurs (1983 hereafter dV83) assuming the Galacto- 
cent ric radius of the Sun to b e Rq = 8.5 ± 0.5 kpc. Sec¬ 


ond, Bahcall & Soneira (1980 herafter B&S) constructed 
a siniilar model that combined a disk and spheroid com¬ 
ponent in order to match observed star counts as a func¬ 
tion of magnitude, latitude, and longitude (rather than 
only the distribution of light across the sky). This work 
yielded global values of Mb = —20.1, My = —20.5 
and B — V = 0.45 (assuming no reddening due to dust 
obscuration; no error estimates were provided). dV&P 
also summarized a series of earlier determinations of the 
Galaxy’s total absolute magnitude made before its mor¬ 
phology and gross stellar structure were well understood. 


— 19.9, or Ml/ 

- -20.2 to -20.5 (Kreiken||1950l 

de Vau- 

couleurs 

1970 

Schmidt-Kaler & Schlosser 

jm 

1 , corre- 


tween 0.3 and 1.0 mag. Additionally, dV83 averaged 
the colors from a set of galaxies believed to have nearly 
the same morphological type as the MW (assumed to be 
Sb/c) in order to infer B — V = 0.53 ± 0.04. 

In the last decade, there has been a growing movement 
to quantify ho w typical the MW is amongst galaxies of 


I 

dl 


its type (e.g., Flynn et al. 2006 Hammer et al. 2007 


Yin et al. 2009D; most ot this work has used the mea- 


work done by|Mutch et al. 

( 

V 

2011 

hereafter Mil) which 

investigates whether the iVI 

V IS located in the so-called 

“green valley” (cf. Mendez et al. 

2011 

Jin et al.| 2014 


gion between the bi modal distribution of red and blue 
galaxie s in the CMD ( Strateva et al.|2001 Blanton et al. 
2003a). After converting van der Kruit’s measurement 
ot Johnson B — V to Sloan Digital Sky Survey (SDSS) 
AB model u — r and placing the MW on a color-stellar 
mass diagram. Mil found that no secure conclusions as 
to the Galaxy’s color could be drawn. 

To help resolve this question, we present in this pa¬ 
per a new method of determining our Galaxy’s global 
photometric properties with dramatically smaller uncer¬ 
tainties. Our technique resembles the “sosies” method 


utilized by Bottinelli et al. 1 

1985 

) and 

de Vaucouleurs & 

Gorwin 

(lyytip. ihe under 

lying 

idea behind that tech- 


brated properties, it can be assumed that they share the 
same luminosity, and hence differences in their apparent 
brightness can be used to determine their relative dis¬ 
tances. Here, we also look for sosies (i.e., analogs) of the 
MW; however, our goal is different, and we take advan¬ 
tage of larger datasets and more sophisticated statistical 
treatments in order to take into account uncertainties 
properly. We derive our results using a method similar 
to that producing the dV83 value; i.e., we average the 
observed properties of galaxies selected as MW analogs, 
though here we carefully account for the systematic bi¬ 
ases that can affect such an approach. 

Essentially, we make the Copernican assumption that 
the MW should not be extraordinary for a galaxy of its 


stellar mass, M*, and star formation rate (SFR), M*. As 
these two properties are very strongly correlated with 
galaxies’ global photometric properties, we first obtain 
a sample of MW analog objects that collectively match 
the stellar mass and SFR of our own Galaxy (taking into 
account the relevant uncertainties). The range of ob¬ 
served photometric properties of galaxies in this sample 
provides tight constraints on our Galaxy’s color and ab¬ 
solute magnitude. With these values determined we are 
then able to accurately determine the MW’s position in 
color-magnitude space. 

Throughout this paper, all SDSS ugriz magnitudes are 
reported on the AB system, whereas all Johnson-Cousins 
UBVRI magnitudes are reported on the Vega system. 
We use a standard ACDM cosmology with Dm = 0.3 and 
Da = 0.7. All absolute magnitudes are derived using a 
Hubble constant of Hq = lOO/i km Mpc“^, therefore 
making them measurements of M — 5 log h. However, in 
order to compare measurements of the Galactic SFR and 
stellar mass, which are measured on an absolute distance 
scale, to those for extragalactic objects measured on the 
cosmic distance scale directl y, we assume a Hubble pa - 
rameter of h = 0.7, following Brinchmann et al. (2004). 


Consequentially, the log M* and log M* values we use for 
external galaxies can be adjusted to different choices of 
Hq by subtracting 21og(h/0.7). For clarity, in what fol¬ 
lows we will explicitly display the h-dependence of all 
quantities we use, as well as explain how our results for 
MW properties change with respect to h. 

We structure this paper as follows. In ^we describe 
our observational data; this includes discussion of our to¬ 


tal stellar mass and SFR estimates for the MW in 12.1 


as well as discussion of the sample of externally measurec 
galaxies we employ in §2.2| In ^ we describe the criteria 
used in order to select the subsamples of SDSS galaxies 
used in this study. In particular, we describ e the selec¬ 
tion of a sample of MW analog galaxies in §3.2[ which 
we use to produce tight constraints on the integrated 
optical-wavelength properties of the Galaxy in |5] In 
we investigate the principal sources of systematic error 
that may arise from our analog-sample selection meth¬ 
ods. We present our final results in 1^ including tables 
of useful photometric properties for the MW. Lastly, we 
summarize this work and discuss its implications in 

2 . OBSERVATIONAL DATA 

In this section, we present a summary of the observa¬ 
tional data we use for this study. We begin by focus¬ 
ing on the total stellar mass, M*, and SFR, M*, of the 
MW. With these parameters in hand, we then describe 
the uniform parent sample of galaxies used, including 
the methods used to measure their stellar masses, SFRs, 
and rest-frame magnitudes. The overarching goal of this 
study is to use this uniformly measured set of galaxies 
to convert our knowledge of the stellar mass and SFR of 
the MW into constraints on its global photometric prop¬ 
erties. The following section will detail how we construct 
a set of MW analog galaxies for that purpose. 

2.1. The Milky Way 


In Licquia & Newman (2015 hereafter LN15), we 
present updated constraints on the total stellar mass and 
SFR of the MW, incorporating the wide variety of mea- 
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surements in the literature. For many of the same rea¬ 
sons that measuring the photometric properties of the 
MW is difficult (cf. 0 , there are a limited number of 
estimates of both Galactic parameters in the literature. 
In order to extract as much information as we can from 
these measurements, which encompass a variety of dif¬ 
ferent observational data and methods, we employed a 
hierarchical Bayesian (HB) analysis method to combine 
all the measurements of a quantity into one aggregate 
result. The HB method allows us to account for the pos¬ 
sibility that any one of the included MW measurements 
is incorrect or has inaccurately estimated errors (e.g., due 
to neglected systematics). The probability of erroneous 
measurements being incorporated into our meta-analysis 
is quantified by the inclusion of hyper-parameters in the 
Bayesian likelihood that characterize the data itself, and 
which we can simultaneously fit for along with the physi¬ 
cal parameters of interest (e.g., M* or M*). The results of 
this study show that the conclusions from an HB analysis 
are robust to many different ways of modeling erroneous 
measurements. 

With the present work in mind, the final results from 
LN15 are normalized so that they can be directly com¬ 
pared to the stellar masses and SFRs of external galax¬ 
ies in the MPA-JHU catalo g. For the SFR of th e MW, 
we capitalize on the work of Chomiuk & Povich (2011), 


which tabulated M* measurements made in the last three 
decades, renormalizing each to a uniform choice of the 
Kroupa broken-power-law initial mass function (IMF; 
Kroupa & WeidnerH2003|) as well as stellar population 
synthesis (SPS) code. Applying the HB analysis method 
to these updated measurements yields a global SFR for 
the MW of M* = 1.65 ± 0.19 Mq yr"!. 

For the total stellar mass of the MW, LN15 apply 
the HB analysis method to nearly 20 independent mea¬ 
surements of the stellar mass of the bulge component 
(including the contribution from the bar) from the lit¬ 
erature, including results from photometric, kinematic, 
and microlensing techniques. For the disk component 
of th e Galaxy, we assume the single-exponential model 
from Bovy & Rix (2013); this is developed from the 
dynarnical analysis of ^16,000 G-type dwarf stars seg¬ 
regated into 43 mono-abundance populations based on 
their position in [a/Fe]-[Fe/H] space, as measured by 
the SDSS/SEGUE spectroscopic survey. Through Monte 
Garlo techniques we are able to simultaneously pro¬ 
duce model-consistent realizations of the bulge and disk 
masses; we sum these two components to yield the to¬ 
tal stellar mass of the Galaxy, M* (the contribution of 
the stellar halo is negligible). The Monte Carlo tech¬ 
niques allow us both to ensure that each bulge mass 
estimate is placed on equal footing and to incorporate 
the current uncertainties in the Galactocentric radius of 
the Sun, Rq. In particula r, we assume the const raints 
of Rq = 8.33 ± 0.35 from Gillessen et al. (2009). We 


show that once the bulge mass estimates are renormal¬ 
ized to the same definition of stellar mass (including 
main-sequence stars and compact remnants, but not 
brown dwarfs), scaled to the same Rq appropriate to the 
measurement technique, and normalized to reflect con¬ 
sistent assumptions about the structure and demograph¬ 
ics of the stellar populations (Kroupa IME and single¬ 
exponential profile disk) then the results from our HB 


analysis are insensitive to models of potential systemat¬ 
ics affecting the data. All of this work culminates in a 
total stellar mass for the MW of M* = 6.08 ± 1.14 x 10^° 
Mq. 

2.2. SDSS Galaxies 
2.2.1. Photometry 

To select a comparison sample of externally measured 
galaxies, we make use of data fr om the Eighth Data 
Release (DR8; |Aihara et ah 2011) of the S loan Digital 
Sky Survey HI (SDSS-Hl; York et al. 2000). DR8 pro¬ 


vides both imaging and spectroscopic data for almost 10® 
galaxies in the local universe, spanning over a third of the 
night sky. Its five broad optical passbands, labeled u, g, 
r, i, and z in order of increasing effective wavelength, 
fully encompass of the GCD-wavelength window. We 
make use of DR8, made available in early 2011, due to 
its best-to-date calibration and reduction of the imag¬ 
ing data. All subsequent data releases from SDSS-HI 
have provided no further refinements for low- 2 ; galaxies 
as studied here. The Photo pipeline processing yields a 
variety of magnitude measurements based on fitting both 
a pure de Vaucouleurs and a pure exponential profile to 
the surface brightness distribution of each object. Those 
quantities labeled “model” reflect the magnitude derived 
from the better of the two model prohle hts in the best- 
measured band (generally r), which is then convolved 
with the object’s point spread function (PSE) in each 
passband to obtain a template for measuring its flux. 
DR8 also provides the magnitude derived from the op¬ 
timal linear combination of the two model profiles that 
best fit the 2D image of any object in each passband, 
again convolved with the object’s PSF; these are labeled 
“composite model magnitudes” or cmodellM The model 
magnitudes are designed to produce the best, unbiased 
estimate of galaxy colors and so we use these to evaluate 
any color properties we discuss below. However, while 
the cmodel magnitudes are not recommended for pro¬ 
ducing galaxy colors, they do reflect the best estimate of 
the “total” flux of a galaxy in each passband. Therefore, 
all absolute magnitudes described below are derived from 
the cmodel measurements. 

We have obtained AT-corrections on all magnitudes in 
the DR8 catalog to rest-frame 2:=0 and 2;=0.1 SDSS 
passbands using the kco rrect v4_2 software package 
(Blanton & Roweis 2007). This entails fitting spectral 
energy distribution (SED) models to the observed ugriz 
extinction- and AB-corrected magnitudes, given the ob¬ 
served redshift, and then using this fit to determine off¬ 
sets between observed quantitie s and magnitudes mea¬ 
sured in rest-frame bands (e.g.,|Hogg et al.||2002|). The 
observed 2 : also provides a luniinosity distance (given 
the cosmology we assume) and hence the distance mod¬ 
ulus, m — M] we use the kcorrect v4_2 software to 
obtain rest-frame absolute magnitudes that are derived 
by subtracting this distance modulus along with the K- 
correction from the extinction- and AB-corrected appar¬ 
ent cmodel magnitude in each band. We obtain galaxy 
colors by taking the difference of two rest-frame abso¬ 
lute magnitudes, but using model magnitudes in place of 

^ See http://www.sdss3.org/drS/algo rithms/magnitudes .p hp 
for further detail, as well as discussions in|Dawson et all] ||2Uld|l. 
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cmodel a s described above. We cho ose to adopt the nota¬ 
tion from Blanton & Roweis (20071 when presenting our 
results: we denote an absolute magnitude for passband 
X as observed at redshift z by ^M^- 

At this point, we also use the kcorrect package to con¬ 
vert each galaxy’s set of SDSS ugriz (AB) magnitudes 
to an equivalent set of Johnson-Cousins UBVRI (Vega) 
magnitudes, as well as their respective AT-corrections. 
This allows us to calculate UBVRI extinction- and 
AT-corrected (cmodel-based) absolute magnitudes and 
(model-based) colors, which we can then analyze in par¬ 
allel to ugriz measurements in order to yield our results 
transformed to the Johnson-Cousin system. As we will 
see in this will be useful for comparing our results 
to the literature, and should be more robust than using 
any transformation equations available that are averaged 
over all galaxy types. 

2.2.2. MPA-JHU Stellar Masses and SFRs 

For a large sample (^10®) of galaxies with spectro¬ 
scopic redshifts from SDSS below 0.7, the MPA-JHU 
galaxy property catalog provides estimates of total stellar 
masses and SFRs. These are currently publicly available 
at http://www.mpa-garching.mpg.de/SDSS/DR7/, and 
are based on SDSS Data Release 7 photometry. However, 
to ensure the greatest possible accuracy in our results, 
for this study we have produced an upgraded version of 
this catalog by recalculating each galaxy’s M* and M* 
using the same algorithms, but applied to the photomet¬ 
ric measurements released in DR8. Hence, our initial 
dataset consists of the subset of galaxies in the MPA- 
JHU catalog that also have photometric measurements 
reduced through the DR8 pipeline. All results presented 
herein are based on our DR8-based M* and M* measure¬ 
ments, which assume a Kroupa IMF. In the following, 
we briefly summarize the Bayesian methodology used to 
produce them. 

Total stella r masses are determined follo wing the same 


philos ophy as Kauffmann et al. (2003) and Gallazzi et al. 
(2005), but by htting models ol SPb to each galaxy’s pho¬ 


tometry instead of using any spectral features. Here, we 
first construct a la rge grid of galaxy models from |Bruzual| 
& Chariot (2003, BC03), encompassing a wide range of 
possible star forrnation histories. Each model produces a 
synthetic spectrum which we convolve with ugriz pass- 
bands to produce model photometry. For each galaxy, we 
then determine the likelihood for each model by calculat¬ 
ing the from differences between fluxes corresponding 
to the model photometry and the observed model magni¬ 
tudes. Adopting flat priors on all model parameters, we 
then calculate the posterior probability for each model 
given t he observatio n s. Th is is most similar to the meth¬ 
ods of Salim et al. (2007), differing in that the latter 
generated sets of input parameters by randomly draw¬ 
ing them from their priors instead of employing a grid. 
Lastly, we integrate our grid of posteriors along all but 
the stellar mass axis in order to produce the marginalized 
posterior PDF for M*. 

SFRs are determined from the technique described 


Brinchmann et al. (2004), but with several improve- 


fluxes measured from their SDSS spectra, after subtract¬ 
ing the continuum and absorption features using the SPS 
spectra from the latest updates to the BC03 libraries. In 
this case, a grid of ^2 x 10® CLOl models are investi¬ 
gated, which make up a four-dimensional grid of metal- 
licities, ionization parameters, total dust attenuations, 
and dust-to-metal ratios. Similarly as described above, 
the resulting grid of posteriors for all models can then 
be integrated over the other three axes to produce the 
marginalized posterior PDF for dust attenuation. This 
is then used to estimate an unattenuated Ha luminosity, 
which is then converted to a SFR using the [Kennicutt] 
(19981 conversion factor. 

'i'his yields a measurement of the SFR of each galaxy 
inside the SDSS 3" hber. To overcome aperture bias, and 
hence produce an estimate of M,^ for the ent i re gala xy, we 


now follow in the footsteps of Salim et al. (2007). This 
requires calculating photometry tor the lignt that falls 
outside of the fiber and fitting stochastic SPS models to 
it; for each galaxy we combine the SFR measured from 
inside and outside of the hber to determine its total M*. 
As a result, the SFRs employ ed herein should ma tch well 
with the “UV” estimates by Salim et al. (2007) for all 


classes of galaxies over the entire dynamical range of M* 
val ues. For a more in-depth discussion, the reader should 


ments. h'or star-forrning g alaxies this entails htting the 


emission line models from Chariot & Longhetti (2001 
CLOl) to their Ha, O II, H/3, U III, i\ II, and S II emission 


see Brinchmann et al. (2004) and the MPA-JHU catalog 
website (listed above). 

Ultimately, our methods produce DR8-based posterior 
PDFs for the log stellar mass and log SFR of each galaxy 
in our SDSS sample. In our discussions to follow, when 
referring to a galaxy’s log M* or log M* we are truly re¬ 
ferring to the mean value measured from the posterior. 
We could equally have used median values for this study, 
as using them instead yields no differences in our results. 
We also calculate the cumulative distribution functions 
(CDFs) measured from each galaxy’s posteriors, and we 
use Vx to denote the value corresponding to the xth per¬ 
centile in the CDF. We then calculate an effective stan¬ 
dard deviation for both quantities as ((P 84 — d’i6)/2- We 
label these as criogM* and criogM., hereafter. We note that 
these effective error estimates are used only to s creen 
galaxies with highly uncertain measurements in §2.2.3| 
and for investigating the impact of Eddington bias on 
our results in §4.1[ and hence are sufficient for our pur¬ 
poses. 

2.2.3. Initial Cuts 

From the sample of galaxies we have described 
so far, we next restrict to a subset of those that 
make up the SDSS main galaxy spectroscop ic sam¬ 
ple (whose overall selection is described in Strauss 
et al. 2002), which includes only objects with good- 
quality, clean measurements. To do so, we take ad¬ 
vantage of the Photo pipeline processing flags and im¬ 
age bitmasks to eliminate problematic objects from 
the full DR8 sample. We first restrict to objects 
that were targeted as main sample galaxies by en¬ 
forcing that the primTarget flag is set to “galaxy.” 
We then reduce to galaxies with good-quality observa¬ 
tions taken from the Legacy target plates by requiring 
the SDSS plate information tags survey, programiName, 
and plateQuality are set to “sdss,” “legacy,” and 
“good,” respectively. We ensure a good-quality detec- 
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tion by requiring that the BINNEDl, BINNED2, or BINNED4 
flag is set for the r-band image. We exclude ob¬ 
jects with any of the following r-band image flag^set: 
SATUR, BRIGHT, BLENDED, NQDEBLEND, DEBLEND_NOPEAK, 
DEBLENDED_TOO_MANY_PEAKS, PEAKCENTER, NOTCHECKED, 
CR, NOPROFILE, MANYPETRO, NOPETRO, PSF_FLUX_INTERP, 
BAD_COUNTS_ERROR,INTERP.CENTER, BAD_MOVING_FIT, or 
DEBLENDED_AT_EDGE. At this point, we also exclude any 
galaxy whose inverse variance (= l/cr^) in absolute g- or 
r-band magnitude is calculated to be 4 mag or smaller 
after AT-corrections, or that has > 1 or criogM* > 

0.5, in order to exclude any object with highly uncertain 
luminosity, color, SFR, or stellar mass (these restrictions 
on property errors exclude only the most extreme outliers 
in the data, comprising ^ 1% of the sample). As a re¬ 
sult of these cuts, the DR8 sample is reduced to 337,331 
galaxies from ^10® with no restrictions applied. 

3. CONSTRUCTING USEFUL SDSS GALAXY 
SAMPLES 

We next trim our set of cleanly measured galaxies from 
the main galaxy sample to produce a uniform subset suit¬ 
able for statistical analyses. In this section we discuss 
the cuts employed to produce two important subsamples 
used in deriving our final results. First, we describe the 
selection of a volume-limited sample, which consists of 
all galaxies lying in a redshift range such that any ob¬ 
ject with both SFR and total stellar mass values similar 
to those of the MW will be included in the SDSS sam¬ 
ple. Next, we discuss our method of identifying a set of 
MW analog galaxies from this volume-limited sample. 
These galaxies, chosen based upon their M* and M* val¬ 
ues, can then be used to estimate the global photometric 
properties of the MW, while the volume-limited sample 
provides the context for discussing the MW’s location in 
color-magnitude space. 


3.1. Selection of a Volume-limited Sample 

The SDSS main galaxy sample ( Strauss et al.||20d2 ) is 
bounded by a limiting Petrosian nmgnitu3e''brr^T7.77 
after correction for Galactic extinction. Of course, this 
introduces a radial selection effect, known as Malmquist 
bias, whereby only the intrinsically brightest galaxies 
are present in the data at large z, whereas less lumi¬ 
nous galaxies at the same redshift will not be targeted 
for spectroscopy. We therefore restrict ourselves to a 
volume-limited sample, ensuring that all galaxies in the 
MW SFR and M* ranges are detected and available for 
selection, regardless of their distances. 

To do this, we first select a sample of MW analog 
galaxies via the process described in the following sec¬ 
tion from the clean main galaxy sample (i.e., ones with 
good-quality measurements, as described in (2.2), but 
with no restriction in redshift. Next, we overlay this 
sample of analogs upon the clean main galaxy sample 
with z > Zniin in ^ — f) VS. color-magnitude space; 

we then increase Zmin until, by eye, the included objects 
fall as faint as the faintest MW analogs, but no more so. 
The corresponding value of Zmin will then serve as the 
maximum redshift for our volume-limited sample. This 


^ See http://www.sdss3.org/dr8/algorlthms/photo_flags. 
php and sources therein tor explanations ot these Hags. 


is true since we expect that the least-luminous galaxies 
in a sample must be at the smallest available redshift 
to be seen at all; or, we can conclude that at any red¬ 
shift below Zmin such faint galaxies would be included 
in a sample, but at any greater redshift they would not. 
Therefore, z„iin corresponds to the upper bound on the 
range of redshifts allowed for a volume-limited sample of 
MW analogs. 

We note that this method is an extension of the stan¬ 
dard procedure generally used to identify volume-limited 


samples of objects (see, e.g., Tago et al. 2010 Tempel 


et al. 2014). Whereas one typically investigates the lu- 


mmosity completeness level as a function of redshift, we 
have extended this analysis to the GMD. In this way, 
we have ensured that all the results we present below 
(i.e., the luminosities and colors we infer for the MW) 
are guarded against Malmquist bias. We also note that 
this process of choosing a z^nin contributes a negligible 
amount of uncertainty to our final results presented be¬ 
low. For example, changing z^in by ±0.005 yields a shift 
in all of our results by < O.OSct. 

Based on investigation of MW analog GMDs for dif¬ 
ferent limiting redshifts, we find that a cut of 0.03 < 
z < 0.09 ensures that all analogs have r < 17.77, so that 
the SDSS magnitude limit has no effect on our results. 
The lower bound on z is used to limit the impact of aper¬ 
ture effects on the properties measured for these galaxies; 
again, a ±0.005 shift in this value yields a < 0.05cr change 
in our results. In addition to applying this redshift cut, 
we simultaneously enforce that all galaxy redshifts are 
measured at high confidence by ensuring that no red¬ 
shift warning flags are set within the SDSS catalog (i.e., 
each has a value z_warning = 0). The resulting volume- 
limited sample includes 124,232 galaxies from the clean 
main galaxy sample. 

In Figure we the show position of the MW with 


la constraints as determined by LN15 in M*-M* space, 
overlaid upon log-spaced contours depicting the density 
of the volume-limited sample. We have highlighted the 
approximate locations of the “main sequence” of star¬ 
forming SDSS galaxies in blue and the region of qui¬ 
escent galaxies in red. We also do this in our color- 
magnitude plots below; note that the relative positions 
of the regions corresponding to these two populations are 
flipped in magnitude space. In the following section, we 
exp lain how we apply the Galactic constraints discussed 
in §2.1| to the volume-limited sample in order to con¬ 
struct a sample of MW analogs, which in total should 
exhibit the same properties as the MW. We can then ex¬ 
amine where these analogs lie in color-magnitude space, 
ultimately converting our knowledge of where the MW 
lies in the M*-M* plane into similar constraints on its 
photometric properties. 


3.2. Identifying Milky Way Analogs 

We now collect a set of galaxies that, as an ensemble, 
can be used to constrain the overall photometric proper¬ 
ties of our Galaxy; i.e., a sample of MW analog galaxies. 
By our definition, these analogs are selected in such a 
way that the distributions of their measured M* and M* 
values match the posterior probability distributions de¬ 
scribing the Galactic M* and M* found in L N15 using 
a HB analysis (these results are detailed in (2.1). To 
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Fig. 1. — Position of the Milky Way in star formation rate (M*) 
vs. total stellar mass (M*) space. Log—spaced contours depict 
the density of a volume-limited sample of SDSS galaxies in the 
range 0.03 < 2 < 0.09. The most likely position (red dot) and 
Icr constraints (purple) for the Milky Way shown here are deter¬ 
mined from a hierarchical Bayesian meta-analysis of the literature 
estimates by LN15. The properties for both the Milky Way and 
the extragalactic sample displayed here reflect consistent assump¬ 
tions about their stellar populations, including a Kroupa IMF, and 
hence should be well guarded from any substantial systematics rel¬ 
ative to one another. We can see the Galaxy is offset from the 
main sequence of star—forming galaxies, hinting that it may be in 
a transitional evolutionary phase. 

do so, we apply a randomized selection procedure to the 
galaxies in the volume-limited sample, as follows. 

We begin by randomly drawing a single value from 
each of the adopted PDFs describing the MW’s M* and 
M*, independently of each other. Ideally, we would like 
to then select a single galaxy from the volume-limited 
sample whose measured properties match these values 
exactly; we could then trivially build a sample of MW 
analogs by repeating this process a large number of times. 
However, in general there will be no galaxies with proper¬ 
ties that match these values perfectly. Therefore, we use 
our pair of values drawn from the Galactic distributions 
as a point of reference in the M*-M* plane and iden¬ 
tify all galaxies from the volume-limited sample that lie 
within a small tolerance window centered on this point. 
We choose this window to be the rectangular box that 
encompasses all values to within x% of the M* and M* 
values drawn. To ensure that the distributions of analog 
properties will still match the fiducial Galactic posteri¬ 
ors, we require that x% is much smaller than the error in 
either of the MW M* and M* results presented in LN15. 
Finally, from the galaxies that lie within our tolerance 
window we randomly select one as a MW analog. We re¬ 
peat this process 5000 times as the first step in building 
our sample, providing us with a set of 5000 MW analog 
galaxies. 

We have employed a tolerance window in our method, 
as opposed to simply selecting the object in the volume- 


limited sample nearest each (M*, M*) pair drawn, to 
maximize the number of unique MW analogs that make 
up our sample. In practice, we find that when using a 
1% tolerance our window encompasses at least one galaxy 
from the volume-limited sample ~75% of the time, and 
typically contains up to eight candidate objects. The 
remaining ^25% of the time we can expand our win¬ 
dow to a 3% tolerance from the drawn M* and M* val¬ 
ues in order to encompass a set of at least one galaxy, 
from which we randomly draw one analog. Given that 
the fractional error in the adopted M,^ and M* for the 
MW is ^19% and ~12%, respectively, we find that the 
exact size of the l%/3% acceptance window is inconse¬ 
quential to this study. We have tested for the impact of 
using broader parameter space window sizes, though still 
small compared to the MW measurement errors, and al¬ 
ways recover the same results (to well within the quoted 
errors). We have also tested for any changes in our re¬ 
sults when selecting analogs by their total specific SFR 
(sSFR) and M* in place of SFR and M*; again, the dif¬ 
ferences are well within the uncertainties. In light of this, 
we have chosen to present the results of using the SFR 
measurements only, as using sSFR introduces substan¬ 
tial covariance (i.e., sSFR correlates strongly with M*), 
whereas the SFR and total stellar mass of the MW are 
determined independently in LNI5. 

Just as it is a problem for observing the MW, dust al¬ 
ters the observed colors and magnitudes of star-forming 
galaxies observed with high inclination angles. We there¬ 
fore exclude objects likely to be edge-on spiral galaxies 
from our Milky Way analog sample (MWAS). Accord¬ 
ingly, from the 5000 galaxies selected initially, we exclude 
all those that have both a best-fit axis ratio b/a < 0.6 
measured from a purely exponential profile fit to the sur¬ 
face brightness density in the i —band, as well as a value 
/deV < 0.5, where /deV effectively denotes the fraction 
of light in the galaxy’s image that is contributed from a 
bulge-like component vs. a disk-like component, again 
generally measured from the 2D r-band image. In ef¬ 
fect, we are excluding all galaxies we have selected that 
are both edge-on and disk-dominated. It is important to 
note that, because we have applied no morphological con¬ 
straints on MW analogs, excluding disks in this way will 
introduce a morphological bias into our sample to some 
extent (i.e., the ratio of bulge-dominated and elliptical 
galaxies to disk-dominated galaxies will increase), an ef¬ 
fect we will need to correct for. Therefore, it is important 
to apply this cut only later in our MW analog selection 
process so that we are able to track the fraction of disk 
galaxies that make it into the sample before and after 
its implementation; knowing these n umb ers will allow us 
to make the proper correction. In §4.2| we will provide 
a more in-depth discussion of this, including how this 
inclination cut was chosen and its impact on our results. 
Ultimately, after removing edge-on disks we are typically 
left with a clean sample of ^3500 galaxies. 

For the particular realization we use in this study, our 
process more precisely yields 3402 galaxies that we will 
use to derive our results below, and which we henceforth 
call the MWAS. We note that, of the analog galaxies se¬ 
lected, only 935 (or ~27%) are unique objects. It is im¬ 
portant to keep duplicate objects so that the distribution 
of property values for the MWAS accurately matches the 
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Parent Sample of Cleanly Measured SDSS Galaxies 

(337,331 objects) 

i 

Discard all objects that do not have all of the 
following properties: 0.03<z<0.09 and 
z_warning = 0 

I 

Volume-Limited Sample 

(124,232 objects) 



•• 


i 

Randomly draw 5,000 pairs of values 
(M*, M*) from the fiducial MW PDFs. For 
each pah', randomly select one object from the 
volume-limited sample with properties within 

1% of these values. If no objects available, then 
retry using a 3% tolerance. 

n 



TI 

ii 

• k 

Discard all objects that do have all of the 
following properties measured from their 2D 
r-band image: f^^y<0.5 and (b/a)^^p<0.6 


r? r.'-i— 



Milky Way Analog Sample 

(3,402 objects) 



Fig. 2.— Flowchart outlining the steps and criteria we use to se lect different samples of SDSS galaxies that we employ in this study. 
This chart summarizes the processes described in sections |2.2H3.2| where more details may be found, including where we obtain or how 
we produce different property measurements for each object. INote th at he re we denote the error in galaxies’ stellar mass and SFR as 
(T(logM*/M 0 ) and crjlogM*/MQyr“^), respectively. As mentioned in [3.2 { 6 /a)exp is the minor-to—major axis ratio obtained from the 
pure exponential profile best fit to an object’s 2D image. 


posterior distributio ns fo r the MW properties we have 
found in LN15 (see (2.1). In practice, we find that if we 


only keep the set of unique objects, the mean M* of our 
sample has a significant bias (~2x the standard error) 
compared to when we eliminate duplicates. The SFR 
distribution is affected less. These biases are avoided al¬ 
together by allowing objects to be selected multiple times 
as a MW analog. 

For convenience, w e show a flowchart in Figure that 


summarizes sections 2.2 |3(2 into a step-by-step proce¬ 
dure that yields all of tne different samples of galaxies 
that we employ in this study. Figure shows the posi¬ 
tions of the MWAS in M*-M* space as red dots, overlaid 
upon the same contours for the volume-limited sample 
as Figure^ The spread of these dots appears broader 
than the Galactic constraints in Figure due to both the 
saturation of color where there are many objects and the 
substantial number of >3cr events to be expected in any 
sample of 3500 numbers; as mentioned above, the size of 
the search box is small in comparison to the spread in 
MW values. 

Figure shows our sample of MW analogs (red dots) 
overlaid on the volume-limited sample (grayscale con¬ 
tours), similar to Figure]^ but now plotted in the '^{g — r) 
vs. ^Mr CMD. Mapping these galaxies from one param¬ 
eter space to the other noticeably increases their scatter 
compared to the underlying population from which they 
were drawn. However, their tight correlation in M*-M* 


space, as expected, produces significant constraints in 
the CMD, providing us with information on what loca¬ 
tions could feasibly be occupied by the MW. We display 
a division line between the red sequence and blue cloud 
regions obtained by taking a line parallel to the slope of 
the red sequence, but offset to the point where contribu¬ 
tions from red sequence and blue cloud galaxies are equal, 
determined based upon the assumption that the spread 
in red sequ ence colors about the center line is Gaussian 
(G. Graves 2012, private communication; cf. Taylor et al. 
2015|) . 'i'he peak density of the MW analogs lies near the 


division line, though many lie far above or below it. 

With the MWAS in hand, we are now ready to cal¬ 
culate new constraints on the photometric properties of 
the MW. For instance, simply calculating the mean and 
standard deviation of the ^{g — r) colors of our sample 
yields ^ 0.72±0.07, and similarly we find ~ —20.75±0.37 
for However, as mentioned above, before present¬ 

ing our final results we must first account for any major 
potential sources of systematic error in our method (e.g., 
the morphological bias introduced from removing any 
edge-on disks from the MWAS), and make the proper 
corrections. We will next discuss these systematics and 
the corrections that they require in the following section. 

We note that for any mean quantity described here¬ 
after, including those provided in our tables of results, 
we are actually using the Hodges- Lehmann (H-L) esti¬ 
mator (Hodges & Lehmann 1963). The H-L estimator 
is a robust measure of the median of the data, which 









































































Fig. 3. — Sample of ~3500 Milky Way analog galaxies (red dots), 
chosen through a random selection process such that they collec¬ 
tively match the same distribution in SFR-M* space as the Milky 
Way (compare to Figures^ andThe grayscale, log-spaced con¬ 
tour lines depict the densuy of a volume-limited sample of SDSS 
galaxies (0.03 < 2 : < 0.09) that encompasses the Milky Way SFR 
and M* ranges throughout this redshift range, but is not affected 
by any limiting magnitude. This is t he sa me sample from which 
the Milky Way analogs are drawn (see ^2.2| for the details). Funda¬ 
mentally, we make the Copernican assumption that the Milky Way 
should not be extraordinary amongst the set of galaxies of simi¬ 
lar stellar mass and star formation rate, and hence some galaxy in 
that set must have closely matching photometric properties; in this 
study we focus on integrated optical-wavelength properties which 
are all but impossible to measure directly. 

is calculated by determining the median value of the 
set {{xi + Xj)/2} for all pairs i,j. For N Gaussian- 
distributed data points with a standard deviation a, the 
mean has standard error while the median has 

uncertainty ^cr/\/0.64fV. The H-L estimator has an er¬ 
ror of ~ct/\/ 0.955-/V, comparable to the mean, but shares 
the robustness to outliers of the median, making it a su¬ 
perior choice in most cases. Hence, using the H-L esti¬ 
mator should reduce the impact of significant outliers in 
our sample, in contrast to the ordinary mean, but it will 
have smaller errors than the median. 

Instead of calculating for all possible pairs, which re¬ 
quires excessive computation time, we bootstrap this es¬ 
timate by choosing a set of random i, j pairs equal to ten 
times our effective sample size (reducing from ~ 6 x 10® 
total pairs down to a much more manageable 34,020 for 
our typical calculations). This introduces a small amount 
of extra uncertainty (= a j \/2 x 0.64 x 34,020) which 
must be added in quadrature to the nominal standard 
error in the H-L estimator. The net result is that our 
estimator yields uncertainties 3% larger than the true H- 
L mean would. This additional uncertainty is negligible 
compared to our overall errors; hence this technique does 
not introduce any measurable amount of potential bias, 
and the bootstrapped H-L estimator in our application 
will still have significantly smaller uncertainty than the 
median value. 



-18 -19 -20 -21 -22 

^Mr — 5 log h 


Fig. 4. — Sample of ~3500 Milky Way analog galaxies (red dots) 
plotted in SDSS — r) vs. ^Mr space. This is the same sample 
of objects shown in Figure i.e., they are selected to produce a 
distribution of star formation rate and total stellar mass values 
matching the probability distribution describing those properties 
for our Galaxy. Again, the grayscale, log—spaced contours depict 
the density of a volume—li mited sample of SDSS galaxies in the 
range 0.03 < z < 0.09 (see j ]3.1^ . For reference, the dashed green 
line shows a simple SP SS color cut dividing the red-sequence and 
blue- cloud regions (G. |Graves| 2012, private communication; see 
i ]3.2| for more details). 

4. SYSTEMATICS 

In this section we discuss the principal systematic er¬ 
rors and biases that could affect the methods applied in 
this study, other than systematic errors in either of MW 
or extragalactic M* and M* measurements, which we de¬ 
fer discussion of to First, we investigate the impact of 
Eddington bias, i.e., the bias resulting from selecting ob¬ 
jects using quantities that are affected by measurement 
errors. We provide details on how we can estimate its 
overall effect, which we then subtract from our final re¬ 
sults. In addition, we analyze the impact of reddening as¬ 
sociated with observing disk galaxies at an inclination on 
the optical properties of our MWAS. We discuss how we 
identify inclined objects in the SDSS measurements, as 
well as our methodology for mitigating reddening or ex¬ 
tinction effects that otherwise, when neglected amongst 
the MWAS, could distort the inferred photometric prop¬ 
erties of our Galaxy. We will demonstrate that, after 
corrections for these effects, the remaining systematic 
uncertainties from these effects are well below statisti¬ 
cal uncertainties. 

4.1. Eddington Bias 

It is important to address how the uncertainties in our 
stellar mass and SFR estimates affect our results. Specif¬ 
ically, we are drawing each MW analog from a small bin 
in SFR-M* space. For any parameter whose intrinsic 
probability distribution function has significant higher 
derivatives (second or beyond), scatter due to errors will 
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Fig. 5. — Modeling the Eddington bias in the Milky Way analog selection method. Left panel: histograms of the mean absolute —band 
magnitudes produced from Monte Carlo (MC) simulations of selecting a new sample of Milky Way analog galaxies as increasing amounts 
of noise are added to galaxies’ total stellar mass (M*) and star formation rate (M*) values. This noise is drawn from a normal distribution, 
with mean of zero and standard deviation determined by the errors in a galaxy’s estimated M* and M*, and is applied p times successively 
before analogs are selected. Since nominal values are all affected by noise, we denote them as the n = 1 case, and so any further degradation 
is marked n = p + 1. Right panel: a least-squares quadratic fit to the four points yielded by subtracting the mean of {^Mr }from 

that of using the distributions from the left panel, as a function of n. {^Mr)i is measured from the mean of our MWAS 

to produce the point at n = 2. We then use this fit to extrapolate the blue datapoint at n = 1, the ordinate of which should reflect the 
difference between the actual measurements, which are affected by Eddington bias, and what would be measured with zero errors, i.e., the 
quantity we desire; this value is subtracted from the observed absolute i—band magnitude of the sample. This same process is applied to 
each absolute magnitude or color considered in this study, and the bias subtracted is listed in Tables Almost always this offset is 

completely subdominant to the statistical errors of our method; the exception is n-band-based color measurements, for which the bias is 
of the same order as, but still smaller than, statistical uncertainties. Even then, the uncertainty in the bias correction is much smaller than 
other sources of error. 


move more objects from bins with more objects to those 
with fewer as opposed to the converse. This causes the 
observed distribution of values with errors to be biased 
compared to the true, underlying distribution. This phe¬ 
nomenon is known as Eddington bias and is very common 
in astronomy; it is the generalized form of the Malmquist 
bias that affects luminosity distributions. For instance, 
since massive galaxies are rare, a galaxy with a large M* 
estimate is more likely to have an actual stellar mass be¬ 
low that value than above, since there are many more 
objects that could up-scatter than down-scatter. As a 
result, in aggregate the M* values of our MW analogs 
should be biased high. Similar effects could affect SFR, 
luminosity, or color estimates. 

To quantify this bias, we consider a statistical exercise 
of perturbing each galaxy’s M* and M* values by Gaus¬ 
sian noise sampled from their estimated errors, and then 
reselecting a set of MW analog galaxies utilizing the per¬ 
turbed measurements. To be specific, we offset the mean 
log M* and log M* values individually for each galaxy by 
a value randomly drawn from a Gaussian distribution 
centered at zero with a standard deviation of that ob¬ 
ject’s tTiogM* or cTj^gM* value. We perform N = 100 real¬ 
izations of this perturbation process, each time selecting 
a new set of MW analog galaxies in the same manner 


described in §3.2[ Galculating the H-L mean for each 
ugriz property for each realization yields a distribution 
characterizing our nominal results with the effects of Ed¬ 
dington bias applied twice, instead of the single impact 
that should affect our standard sample. We bootstrap 
this distribution of doubly biased values t o measure th e 
mean H-L mean and its standard error (Efron 1979). 
We then repeat this exercise but applying the noise 2, 3, 
or 4 times consecutively before selecting a sample. This 
yields distributions of the mean property of interest af¬ 
ter repeatedly applying the bias in our method n times; 
the mean of this distribution we denote /i„. For clar¬ 
ity, note that we consider the actual Eddington bias in 
our standard MWAS as the first (n = 1) application, and 
thus distributions of ugriz properties yielded from p suc¬ 
cessive perturbations of the M* and M* values by their 
errors are labeled n = p + 1 in our plots and discussion 
below. Figure displays examples of this analysis for 
°Mr. 

To estimate the Eddington bias in each property we 
then plot the difference between the means of the n and 
n — 1 values of a given parameter as a function of n; i.e., 
Mn — Rn-i versus n. We then perform a least-squares 
quadratic fit to these four data points, incorporating the 
error estimates from our bootstrap analysis. We use the 
resulting curve to extrapolate to n = 1, whose ordinate 
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corresponds to the offset in the “mean” of a given prop¬ 
erty between when Eddington bias affects our sample of 
MW analog galaxies and when it does not. This value is 
then subtracted from the observed mean for that prop¬ 
erty of the MWAS. In Figure|^we show what the results 
of this exercise typically look like, again adopting as 
an example. 

In order to calculate the uncertainty in our estimate 
of the Eddington bias, we construct the covariance ma¬ 
trix for the coefficients of a least-squares quadratic fit, 
A + Bn + Cn^. We are interested in the cr of the point at 
n = 1; this simply reduces to the square root of the sum 
of all elements of the covariance matrix. We note that if 
the uncertainties in our stellar mass estimates were pri¬ 
marily due to photometric errors, this treatment would 
be incorrect, because if an object had (say) a higher- 
than-actual estimated M* value, it would also have a 
too-bright Mj.. However, this does not appear to be the 
case; we find that stellar mass errors are > 5x larger 
than would be expected from SDSS photometric errors, 
so other sources of uncertainty clearly dominate, and we 
can safely treat absolute magnitudes and stellar masses 
as statistically independent. 


4.2. Inclination Reddening 


As mentioned in §3.2[ we have removed any edge-on 
disk galaxies that originally were included in the MWAS. 
This is because, just as reddening and extinction affect 
observations though the disk of the MW, they also al¬ 
ter measurements of external spiral galaxies with their 
disks aligned along our line of sight. The SED measured 
for such an edge-on disk galaxy will be significantly dis¬ 
torted; we will receive a much smaller fraction of blue 
light than when observing face-on, while the flux of red¬ 
der optical light detected will be much less affected. This 
means that the more edge-on a disk galaxy in the SDSS 
sample is, the less representative our observations will 
be of its intrinsic photometric properties (see, e.g., Un- 
terborn &: Ryden |2008[ Mailer et al.||2009 Saliin et IlT 
2009|). It is important, therefore, to ensure the proper- 
ties of the MWAS are not skewed by this effect. In this 
section we discuss quantitatively our method of choosing 
the appropriate inclination threshold for the MWAS, as 
well as the unwanted side effect that it creates, namely 
morphological bias, which we will also need to correct 
for. 

DR8 provides measurements of each galaxy’s ratio of 
semiminor to semimajor axis, 6/a, as determined from 
the exponential profile best fit to its 2D image (labeled 
abExp in the DR8 catalog). Low values of h/a indicate 
that the galaxy’s image has high eccentricity. Addition¬ 
ally, when calculating cmodel magnitudes (see ^2.2[ ), the 
weight of the de Vaucouleurs profile in the best-ht lin¬ 
ear combination with the exponential profile matched to 
the object’s image is recorded as /deV (or alternatively 
fracDeV). Essentially, this quantifies the fraction of the 
total light in the 2D image of the galaxy that is well fit by 
a Sersic index of 4 as opposed to 1. For our purposes, we 
classify any object with /dev > 0.5 as a bulge-dominated 
or elliptical galaxy, and any with /deV < 0.5 as a disk- 
dominated galaxy. Our objective is then to use these 
two parameters to effectively eliminate edge-on disk- 
dominated objects (i.e., ones well fit by an exponential 
profile and that appear to have a high inclination angle) 


from being selected as part of the MWAS to prevent any 
systematic offsets in our results due to reddening. We are 
not concerned about including early-type galaxies with 
low 6/a in the sample, as they contain little cold gas and 
dust, so extinction effects are comparatively minor for 
them. 

We have tested the impact of this cut by varying the 
minimum allowed axis ratio for disk-dominated galaxies 
in our sample and measuring how the mean ^{g — r) color 
is affected. To do so, we again employ a process of select¬ 
ing a set of 5000 MW analog galaxies identical to that 
explained in §3.2[ but unlike before, we do not yet apply 
any constraints on the 6/a or /deV values in the sample. 
In Figure we display how the integrated color of this 
set changes as we impose constraints on 6/a and /deV in 
three different ways. 

First, the upper (red dashed) curve shows the effect 
of removing only disk galaxies (i.e., ones having /dev < 
0.5) as we systematically increase the minimum allowed 
axis ratio. Initially, as we increasingly remove the lowest 
6/a (most inclined) disks we see the mean — r) color 
shifts blueward, as expected. However, once we increase 
our threshold to remove disks with b/a 0.35, we find 
an unwanted side effect: increasingly removing the disk 
population gives increasing weight to bulge-dominated 
and elliptical populations, yielding a trend toward net 
redder color. 

Second, we investigate a scenario that attempts to 
avoid this problem. The lower (dashed-dotted blue) 
curve shows the effect of removing any object regard¬ 
less of its type (i.e., its /deV value) as we systematically 
increase the minimum allowed axis ratio. We note that 
at 6/a > 0.35 bulge-dominated and elliptical galaxies 
outnumber disk-dominated galaxies at a ratio of ~3:2; 
this is the case for the entire volume-limited sample, as 
well as MW analogs. Hence, we find that in this case 
the mean color of the sample becomes increasingly bluer 
as we increase our minimum allowed 6/a threshold above 
0.35. However, since this trend does not stabilize as we 
push our threshold higher, it is likely that we are in¬ 
creasingly oversampling disk-dominated objects due to 
preferentially discarding the more prevalent elliptical and 
bulge-dominated galaxies in this regime. 

Lastly, we present an alternative treatment shown by 
the middle (solid black) curve. Here, we remove the same 
disk galaxies that we do in the first scenario, leaving all 
bulge-dominated and elliptical types initially selected in 
the sample. However, we now calculate a weighted mean 
quantity, where we ensure that the contribution of disk 
galaxies remaining, after applying any minimum thresh¬ 
old on 6/a, is equal to that from the disks present in the 
sample before any cut. In other words, we calculate the 
mean property of our filtered sample of galaxies by mul¬ 
tiplying the contribution from the remaining disk types 
after selection by a weighting factor W = A^^si°s'^®/A^dfs£j 
where and A^^fsks represent the number of disk 

galaxies in the sample before and after applying this cut, 
respectively. For instance, for ^{g—r) color our estimator 

















11 



Fig. 6 .— Mean ^{g—v) color of our Milky Way analog sample as a 
function of the minimum allowed axis ratio of the objects included 
in three different scenarios. The red dashed curve is the result of 
removing only disk galaxies (i.e., objects having /deV < 0-5) based 
on the axis ratio cut given on the a:—axis. The blue dashed-dotted 
curve is the result of removing any galaxy, regardless of type, based 
on the minimum ratio cut. The solid black line reflects the result 
of removing only disk galaxies again, but in this scenario we add 
extra weight to the contribution from disk galaxies remaining af¬ 
ter the cut to the overall mean ^{g — r) color to correct for the 
objects removed. Excluding only disk galaxies initially causes a 
trend toward net bluer color as edge-on systems are removed, but 
eventually this trend reverses toward redder average color due to 
increasingly oversampling the elliptical population. However, in 
the second case, since spheroids outnumber disks at h/a > 0.35, 
throwing away any galaxy above this minimum allowed threshold 
means discarding more spheroids than disks, giving extra weight 
to the blue population. Therefore, we chose the last scenario (a 
weighted mean) as our fiducial method, as it provides a stable be¬ 
havior over a large range of reasonable choices of cutoff for 5/a, 
and hence appears robust to such morphological bias. That is, 
the slope of the black curve is shallow enough that moving from a 
minimum allowed b/a of 0.4-0.8 would cause a < 0.01 magnitude 
change; b/a > 0.6 is our fiducial cut. The other scenarios pro¬ 
vide a much more crude and extreme way of removing inclination 
reddening from our sample; we note that the offset of the other 
curves from the black curve at a: = 0.6 is still < 0.02 mag, which is 
subdominant to the statistical error (~0.06 mag). The analogous 
plot for ®r-band absolute magnitude yields a similar conclusion, 
and so we adopt the weighted-mean scenario as standard for all 
quantities. 


reduces to 
(°(5- d) = 

(°(ff ~ ^)diskDH-L + -^ellipticals (°(g ~ ^')ellipticals)H-L 
+ -^ellipticals 

( 1 ) 

where -Veiiipticais is the number of ellipticals in the 
MWAS, (°(g — »')diskDH-L is the mean — r) color of 
disk galaxies after the h/a and /dev cuts are applied, and 
(°(g-r)eiiipticais)H-L is the mean °{g-r) color of ellipti¬ 
cals in the MWAS. The H-L subscript here denotes that 
we are truly using the Hodges-Lehmann estimator of the 
mean. 


By comparing the different curves in Figurej^ it is clear 
that the weighted mean is favorable over the other sce¬ 
narios for two reasons. First, the slope of the weighted- 
mean curve is shallow enough that moving from a cut of 
b/a > 0.4 to b/a > 0.8 makes a < 0.01 magnitude dif¬ 
ference in integrated *^(5 — c) color. Hence, compared to 
the other scenarios, this prescription results in more sta¬ 
ble values of mean °(g — r) over basically all reasonable 
threshold values of b/a. Second, the position of this curve 
is between the other two, indicating that we are avoid¬ 
ing giving too much weight to either of the red or blue 
populations, and hence limiting the impact of any mor¬ 
phological bias. We note that simply cutting out whole 
classes of galaxies based on their axis ratios, as in the 
first two scenarios, provides a more extreme way of deal¬ 
ing with inclination reddening. The difference measured, 
however, between the curve for our weighted-mean sce¬ 
nario and either of the other two scenarios is < 0.02 mag; 
all of these differences are significantly less than the sta¬ 
tistical errors in °(g —r) color of the sample (~0.06 mag). 
Based on the results of this exercise, we eliminate any of 
the original 5000 galaxies selected for the MWAS hav¬ 
ing both b/a < 0.6 and /dev < 0.5, and we report the 
reweighted mean property to minimize the impact of any 
reddening of the sample due to inclination. As noted 
above, this cut typically removes ^30% of a MWAS re¬ 
alization. 

We note that if we were to avoid edge-on disks entirely 
during our selection of MW analogs (they have been se¬ 
lected based only on their M* and M* values), then our 
results would suffer from the same morphological bias 
that is evident from the red dashed curve in Figure 
but with no way of knowing to what quantitative extent. 
Therefore it is important to make a correction only at this 
stage in the analysis — i.e., calculating our results via a 
weighted mean after removing edge-on disks — so that 
we may correct for both systematic effects, namely incli¬ 
nation reddening and morphological bias. Fortunately, 
we have found that the variations due to such effects are 
far below the random uncertainties. 

One other option would be to correct for inclination- 
related reddening on an object-by-object basis, rather 
than trimming an d reweighting the sa mple as was done 
here. For instance, Mailer et al. (20091 provide formulae 
for converting inclihation-dependent observed quantities 
into intrinsic ones. When applied to the color-magnitude 
relationship their methods transform the SDSS blue-to- 
red galaxy ratio from 1:1 to 2:1 in the absolute magni¬ 
tude range of —22.75 < Mk < —17.75. Converting to 
intrinsic properties could eliminate the need for inclina¬ 
tion cuts altogether and allow us to utilize a larger subset 
of the SDSS main galaxy sample, though at the cost of 
adopting a particular model for extinction corrections. 
In any case, we point out that the difference between our 
choice of correction and the extreme limits displayed in 
Figure]^ is a ^0.02 mag shift, which turns out to be a 
factor of 3 times smaller than the uncertainties in our 
final results. In addition, the slope of this curve becomes 
very shallow beyond a minimum axis ratio of 0.4, so any 
reasonably chosen cut would yield negligible change to 
our results. Similar analyses have demonstrated that 
this same method works well for all colors considered as 
well as for correcting extinction in absolute magnitudes. 
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Fig. 7. — Distribution of stellar masses, M* {left panel), and star formation rates (SFRs), M* {right panel), for our Milky Way analog 
sample at different stages of our analysis procedure overlaid upon the Galactic posterior probability distribution functions (blue solid 
curves) used for selecting them. The red dashed line shows the original sample of 5000 analogs drawn before any cuts are applied. The 
black dashed-dotted curve shows the remainder of the original sample after removing those that appear to be edge-on disk-dominated 
systems, whose inclusion would otherwise systematically redden our results for photometric properties, and then renormalizing to reflect 
the original sample size of 5000. The green dash—triple—dotted curve shows those objects that make up the black dashed-dotted curve, but 
reweighted to correct for any morphological bias (i.e., the oversampling of bulge-dominated vs. disk-dominated objects; see Figure]^ that 
our cuts introduce; see Equation Jljl. 


Overall, we expect that any alternative prescription for 
inclination would have inconsequential impact on the re¬ 
sults of this study. 

Lastly, in Figurej^we show the distribution of M* and 
M* values for MW analogs compared to the posterior dis¬ 
tributions used for selecting them. This includes distri¬ 
butions for the original sample of 5000 before any cuts, 
the 3402 galaxies remaining after removing those with 
b/a < 0.6 and /dev < 0.5, and the reweighted distribu¬ 
tion of those 3402 objects in congruence with Equation 
([^. Where necessary, we have renormalized each distri¬ 
bution to reflect a total sample size of 5000 objects. In all 
cases, we hnd that the mean and standard deviation of 
our sample match those of the posterior distribution and 
have confirmed tha t they are Gaussian-distribu ted via a 
Q-Q plot analysis ( Wilk fc Gnanadesikan|196'8 ). Hence, 
we find that our treatment of inclination reddening and 
morphological bias does not compromise the fundamen¬ 
tal design of our MWAS in M*-M* space. 


5. RESULTS 

With the MWAS assembled and major systematic er¬ 
rors accounted for, we are now able to produce a com¬ 
prehensive outside-in portrait of our Galaxy. Table 
presents the inferred photometric properties we deter¬ 
mine for the MW in rest-frame z=0 SDSS passbands, 
and likewise Table presents rest-frame z=Q.l SDSS 
passband results. Ine values shown are calculated as 
the weighted ( Hodg es-Lehmann estimator of the) mean 
as described in §4.2 and have been corrected for Edding¬ 


ton bias as detailed in 14.1 Each row is calculated in¬ 


dependently of any other table entry; for instance, we 
utilize the full distribution of — r) amongst the MW 
analogs, rather than deriving this value by subtracting 
^Mr from °Mg (this is also due to our colors being de¬ 
rived from model magnitudes, whereas absolute magni¬ 
tudes are based upon cmodel). For reference we list the 
inherent Eddington bias that has been subtracted in jux¬ 
taposition to each corrected value. 

In addition, we tabulate the derivative of each property 
with respect to total stellar mass and SFR. This is ac¬ 
complished by offsetting the distributions we assume for 
the Galactic M* and M* by ±0.1 times their respective 
errors and redoing our analyses. Along with our fidu¬ 
cial results, this provides three data points to which we 
fit a quadratic Lagrangian-interpolation polynomial, and 
then calculate its derivative at the central datapoint. We 
choose an offset of ±0.1cr so that the resulting Galactic 
range in M*-M* space does not require selecting a new 
volume-limited sample of objects. 

As discussed in |l] all extragalactic measurements of 
logM* and logM*, measured on the cosmic distance 
scale, can be converted to reflect different values of the 
Hubble constant by subtracting from them 2 log(/i/0.7), 
effectively shifting them relative to the MW’s position 
in this parameter space. If we were to instead add this 
quantity to the Galactic logM* and logM* values, the 
change in our results would be identical; this allows one 
to calculate how the absolute magnitudes and colors we 
calculate for the MW change for different h using quan- 
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TABLE 1 

Photometric Properties for the Milky Way: Rest-frame z=0 SDSS Passbands 


Property 

Corrected Value 

Bias Removed 




(mag) 

(mag) 

(10-1° mag Mq^) 

(mag yr) 

— 5 log h 

iy.lD_Q 4y 

0.240 ±0.014 

-0.05 

-0.32 

^Mg — 5 log h 

-20.36t°'^( 

0.142 ±0.011 

-0.11 

-0.40 

^Mr — 5 log h 

-21.00+gj8 

0.134 ±0.009 

-0.11 

-0.48 

^Mi — 5 log h 

01 oy+O.SS 
_0.36 

0.120 ±0.009 

-0.14 

-0.49 

— 5 log h 

-21.56+°j® 

0.126 ±0.009 

-0.15 

-0.39 

1 

0 

2.0431°;;®® 

0.090 ± 0.0060 

0.07 

-0.02 

°{u-- g) 

i-358l°;J°| 

0.077 ± 0.0047 

0.06 

0.02 

°{g-r) 

0.6821°;°®® 

0.015 ±0.0017 

0.03 

0.00 

°(r — «) 

U.^yD_Q Q4g 

0.012 ±0.0012 

0.01 

0.00 

O(i-z) 

f) 291 +0 043 
u.^yi_Q 041 

-0.001 ± 0.0009 

0.01 

-0.04 


Note. — The Eddington bias estimated for each band , as described in j ]4.1| is listed in Column 3. This is subtracted from the mean 
property measured from the MWAS, as discussed in j4.2| (see Equation 0): m order to produce the corrected value listed in Column 2. 


titles given in Tables For example, 

d(°Mr - 5 log h) d{°Mr - 5 log h) dM* 


dh 


aM* dh 
d(°Mr - 5 log/i) dM* 




dh 


( 2 ) 


To be explicit, this means that calculating absolute 
magnitudes using a different value of h (where h = 
7ro/(100km s“^Mpc“^) has been used) would shift both 
the positions of the MW and the volume-limited sample 
together in unison along the absolute-magnitude axis of 
any CMD we show. However, calculating M* and M* val¬ 
ues using a different value of h (where h = 0.7 has been 
used) would shift the position of the volume-limited sam¬ 
ple relative to the MW’s position in the CMD; the size 
of this effect can be estimated using Equation ([^. For 
instance, for the MW M* and M* values we have used 
along with the values in Table we find that a ±0.05 
shift in h corresponds to a ~ ±0.05 magnitude shift in 
^Mr — 5 log h and no shift in ° (5 — r). Therefore, we would 
expect that any reasonable difference between the true 
value of h and 0.7 will yield negligible changes in our 
results (well below the measurement uncertainties) and 
the conclusions we draw from them. 

Similarly, the results for the total stellar mass of the 
MW found in LN15 would be changed if any adjustments 
are made to the absolute distance scale (see Table 6 of 
that paper). Predominantly this manifests in changes 
to the Galactocentric radius of the Sun, i?o; LN15 con- 
servatively used 8.33 ± 0.35 kpc based on the work of 


Gillessen et al. (2009). Firstly, since we found that 
dM^/dRo = 3.09 x lO^o M© kpc"^ for the MW, the 
impact of a change in Rq can be obtained by replacing h 
with Rq in Equation ([^ (note ^M^,/^Ro = 0). 

Uncertainties in i?o oominate the error budget in our 
M* model. We find tha t if we were to instead adop t Rq = 

, yield- 


(2015 


len the total 


8.36±0.11 kpc based on Chatzopoulos et al. 
ing a ^69% decrease in uncertainty in Ko; 
stellar mass from LN15 becomes M* = 6.18 ±0.50 x 10^° 
Mq, corresponding to a net ^57% decrease in M* un¬ 
certainty. This ultimately yields a ^20% decrease in Mr 


uncertainty, while also causing the MW analogs to lie 
along a tighter trend in color-magnitude space. Ulti¬ 
mately, as our knowledge of the structure of our Galaxy 
improves (e.g., by measurements from Gaia), our meth¬ 
ods should be able to more strongly constrain the MWs 
location in CMDs. In contrast, the current uncertain¬ 
ties in the Galactic SFR have a negligible effect on the 
constraints on the photometric properties derived in this 
paper. This is because the uncertainty in the MW stel¬ 
lar mass is a significant fraction of the range of stellar 
masses amongst galaxies of comparable SFR, while the 
SFR uncertainty is a ~7x smaller fraction of the range 
of SFRs at fixed mass. We note also that the evolution of 
galaxies since z ~ 0.1 appears to have negligible effect on 
our results; e.g., limiting our analysis instead to objects 
at 0.045 < z < 0.075 yields differences in our results that 
are much smaller than the errors. 

Figure]^ now shows the position of the MW corrected 
for Eddington and inclination bias, as listed in Table 
as a red dot in rest-frame SDSS °(g — r) vs. °Mr 
space; it is overlaid upon log-spaced density contours 
for the volume-limited sample. The purple ellipse dis¬ 
plays our Icr confidence region, accounting for the covari¬ 
ance between color and absolute magnitude; this yields 
a vast improvement in constraining how our Galaxy fits 
amongst the extragalactic populat ion compared to the 


previously best Icr constraints from van der Kruit (1986 
gray dashed-dotted lines). For convenience, we have 
highlighted the red-sequence and blue-cloud regions of 
this diagram (flipped in position compared to M*-M* 
space, since higher-SFR galaxies are bluer). We see that 
MW’s position straddles the boundary between these 
two populations, with a chance that it lies in the core 
of the red sequence or redder. In addition, given that 


galaxv population 

(IStrateva et al. 

2001 

2003a Wong et al. 

[2012 S14), we 

see tn 


Blanton et al. 


- 0.37 


the brightest spiral galaxies in the local universe, while 
its integrated color of ^(5 — r) = 0 . 682 to;o 56 ranks it 
amongst the reddest as well. 

Lastly, we produce an updated plot equivalent to Fig- 
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TABLE 2 

Photometric Properties eor the Milky Way: Rest-frame z=0.1 SDSS Passbands 


Property 

Corrected Value 

Bias Removed 

a/9M* 



(mag) 

(mag) 

(10“^® mag Mg^) 

(mag yr) 

o -51og/x 

-18.8SlO;« 

0.271 ± 0.012 

0.02 

-0.38 

° - 5 log 

-20.07l°t! 

0.168 ±0.011 

-0.10 

-0.4S 

O-^Mr - Slogh 

-20.78lOj^ 

0.130 ± 0.009 

-0.10 

-0.42 

o iMi - Slogh 

— 21 16+0-38 

zi.iO_o_37 

0.134 ± 0.009 

-0.13 

-0.S9 

- Slog 

—21 41+0-39 

0.124 ± 0.009 

-0.14 

-0.3S 

o 

1 

Q oni 

Z.ZUi_Q 172 

O.IOS ± 0.0072 

0.10 

0.08 

o 

1 

1 4iq+0-l24 

0.074 ± 0.00S2 

0.07 

0.04 

0.1(g _ 

n yQo+O.OSl 

U. / o^_0,063 

0.031 ± 0.0027 

0.02 

-0.03 

°'^(r — i) 

0 ‘iqn+0-046 

U.oyU_Q Q42 

0.001 ± 0.0010 

0.01 

-0.01 


n 275+0-047 

U.Z/0_Q Q4g 

0.010 ±0.0011 

0.01 

-O.OS 


Note. — The Eddington bias estimated for each band , as described in j ]4.1| is listed in Column 3. This is subtracted from the mean 
property measured from the MWAS, as discussed in j4.2| (see Equation 0): m order to produce the corrected value listed in Column 2. 



Fig. 8.— Eddington-bias—corrected position of the Milky Way 
in SDSS ^{g — t) vs . ^Mr color—magnitude space (red point and 
purple la ellipse). For comparison, we show in gray dashed-dotted 
lines the previously best la constraint directly measured by vdK86, 
converted to th e SPSS AB magnit ude system via transformation 
equations from |Cook et al.| (|2014||. In order to place this mea¬ 
surement in this plot we subtract 5 log(/i/0.7), allowing it to be 
directly compared to the SDSS sample. Log-spaced contours show 
the density of galaxies in our volume-limited sample; we shade the 
core of the red-sequence and blue-cloud regions in red and blue, re¬ 
spectively, and show the same green dashed color division line as in 
Figure]^ Until now, the Milky Way’s position has remained highly 
uncertain in this parameter space. Our new measurement dramati¬ 
cally improves our knowledge of how the Galaxy compares to others 
in the local Universe; we likely straddle the division between the 
blue-cloud and red-sequence populations, or the so-called “green- 
valley” region of this diagram. This ranks the Milky Way amongst 
the brightest and reddest spiral galaxies still producing new stars 
today. It may well be in a transitional evolutionary phase where 
star formation is dying out. 


ure 1 of Mil by showing our constraints on the MW’s 
position in ^(u — r) vs. M* space, where the green valley 
becomes stretched out and more distinguishable. Here, 
we have highlighted the green-valley region based on two 
different definitions. First, the dark green region shows 
the division line empiri cally derived for SDSS galaxies 
by Baldry et al. (2006|) with an offset of ±0.1 mag in 
the vertical direction (the definition of the green valley 
employed by Mil), which matches well with the den¬ 
sity contours for our volume-limited sample. Second, 
the light green region shows a definition based upon cor- 
recting all SDSS galaxie s for dust effects, as defined by 


Schawinski et al. (2014 hereafter S14). In the second 
case, many of the intermediate-color objects are blue 
galaxies that are both dusty and viewed edge-on, and 
so switching to intrinsic (face-on) properties moves this 
population blueward in the plot, effectively thinning out 
and expanding the green-valley region more. Given that 
our measurement of the MW’s position in this space is 
effectively face-on, the green-valley definition from S14 
provides a suitable comparison. 

Compared to the prior constraints (gray dashed-dotted 
lines), we are in a much better position to now iden¬ 
tify where the MW lies relative to other galaxies. In 
particular, our Galaxy appears to be entering, if not al¬ 
ready a part of, the green-valley region where objects are 
expected to be in a transitional phase; here star forma¬ 
tion is quenching by some mechanism(s); consequentially, 
green-valley galaxies are moving on a traject ory toward 


the red sequence (for more detail see, e.g., jGongalves 
eraI]|2012||Fang et d.||2012| S14). 


i'ablesqS] and [Twe present our results for MW 
3rties transform^ 1 


in 

properties transform^ to the Johnson-Cousins passband 
system. As a reminder, these values have been calcu¬ 
lated in an identical fashion to those listed in Tables 0 
and 1^ but after transforming each MW analog’s set of 
SDSb ugriz magnitndes to an equivalent set of Johnson- 
Cousins UBV RI measurements using the kcorrect soft¬ 
ware package (Blanton & Roweis 20071. This entails 
calculating magnitudes in UBVHI passbands from the 
linear combination of template galaxy SEDs from BC03 
that best fits the observed SDSS ugriz photometry on 
an object-by-object basis, and hence should provide 
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TABLE 3 

Photometric Properties for the Milky Way: Rest-frame 2:=0 Johnson-Cousins Passbands 


Property 

Corrected Value 

Bias Removed 

a/9M* 




(mag) 

(mag) 

( 19-10 Mg^) 


(mag Mgl yr) 

— 5 log h 

-20.02tol° 

0.232 ±0.011 

- 0.10 


-0.61 

s — 5 log h 

^U.U/_o 44 

0.173 ±0.009 

- 0.12 


-0.49 

^My — 5 log h 

—20 74+^-37 

ZU./‘4:_o 39 

0.132 ±0.008 

- 0.12 


-0.37 

^Mji — 5 log h 

-21.26/g:4° 

0.131 ± 0.007 

- 0.11 


-0.16 

^Mj — 5 log h 

-21.84/°;36 

0.125 ±0.008 

- 0.12 


-0.34 

°{U-V) 

0.89010;)48 

0.094 ± 0.0055 

0.02 


-0.14 

°{U-B) 

u.i^y.Q 070 

0.063 ± 0.0038 

0.01 


-0.04 

°{B-V) 

n 744+®-^^^ 

'^^-0.054 

0.028 ± 0.0022 

0.01 


-0.05 

°{V-R) 

n 541+0.046 

U.04i_Q Q42 

0.005 ± 0.0008 

0.00 


- 0.02 

°{R-I) 

u.oyo_o 049 

0.007 ± 0.0009 

0.01 


0.01 


Note. — Values in this table are determined from analyzing the distributions of properties for Milky Way analogs, but after transforming 
SDSS ugriz measurements to Johnson-Cousins UBVRI-equivalent values on an object—by—object basis using the kcorrect software. As 
a reminder, UBVRI magnitudes are on the Vega system, whereas ugriz magnitudes are on the AB system. 


the most accurate transformations. A viable alternative 
would be to apply the empirica l color transformations 
provided by |Cook et al.| (|2014|) directly to our results 
in Tables [^and though these equations represent the 
mean traiisformations between the two passband systems 
averaged over galaxies with a range of morphologies. 


SFRs, etc. Nevertheless, we find that applying the [Cook 
|et al.|20T4| transformation equations to our SDSS results 
produces estimates on the Johnson-Cousins system that 
are quite similar to our nominal values determined using 
kcorrect. The differences are almost always at the 0.1- 
0 . 3(7 level (including for mass-to-light ratios, which we 
discuss next), the one exception being °{U — B), where 
the two meth ods agree at the 0.75 ct level. Note tha t we 
have used thelCook et al.| equations to transform the 
der Kruit ( 1986| ) result in Figure 


van 


In addition to the photometric properties presented in 
Tables we also provide in Table new estimates 
of globm stellar mass-to-light ratios, T*, for the MW 
for all SDSS and Johnson-Cousins passbands in the z=0 
and z=0.1 rest frames. These are calculated from the 
full distribution of T* values for the MWAS, in the same 
manner as we calculate photometric properties. To do so, 
we first calculate the stellar mass-to-light ratio for each 
MW analog in passband x in the rest frame of redshift 2 
as 


= M* X loO-4(CM-+51og(O.7//i)-"M,,,0) ^- 1 ^ (- 3 ) 

where and Lq represent the absolute magnitude 

and luminosity, respectively, of the Sun, which we cal¬ 
culate using the k_solar_magnitudes routine from the 
kcorrect package. We note that Equation (1^ is written 
to make it clear that we have converted absolute magni¬ 
tudes to reflect h — 0.7 and hence be on the same scale 
as our M* values; however, it should be noted that T* is 
intrinsically a cosmology-independent quantity. For in¬ 
stance, if we now chose to rescale quantities from h = 0.7 
to 0.8, the right hand side of Equation (§ would gain a 
factor of (0.7/0. 8 )^ for the change in stellar mass and 
a factor of 10 °- 4 ( 5 iog(o. 8 /o. 7 )) j-pj. change in luminos¬ 
ity, which cancel. Next, we use Equation Q, replacing 


— r) with ^T*, in order to produce our weighted- 
mean estimate. Lastly, we multiply this by a factor of 
10-0-4®^ where B is the Eddington bias correction, which 
is listed in Tables [l}]4 In this way, our results are cor¬ 
rected for Eddingtonbias, inclination effects, and mor¬ 
phological bias (all subdominant to the random errors), 
consistent with all other properties presented here. 


6 . SUMMARY AND DISCUSSION 

This paper has focused on determining the global pho¬ 
tometric properties of the MW to facilitate comparisons 
to observations of other galaxies. In LNI5 we have de¬ 
rived a new, highly constrained stellar mass and SFR for 
the MW using a multitude of independent results from 
the literature, which encompass many different methods 
and datasets. We then identified a set of SDSS galaxies 
analogous to the MW, whose distribution of SER and 
total stellar mass values match the probability distri¬ 
butions for these quantities (given uncertainties) of the 
MW. These two quantities are strongly correlated with 
a gala xy’s luminosity and color (see, e.g.. Bell & de Jong 
20011, so a galaxy that matches our Galaxy in stellar 


mass and SFR would also be expected to have a similar 
overall SED. We then determine the range of photometric 
properties of these galaxies, allowing us to constrain MW 
properties in a manner that is largely robust to the ef¬ 
fects of Galactic extinction (unless the MW is so unusual 
that it has no true analogs amongst the set of galaxies 
matching its M* and M*). We have accounted for the 
Eddington bias involved with selecting galaxies based on 
their SER and stellar mass, and tested the impact of red¬ 
dening effects on this sample. In we have provided a 
full tabulation of useful MW photometric properties. 


6.1. Comparisons to Earlier Color Measurements 

Overall, the results from our MW analog-based analy¬ 
sis method compare well with literature estimates of the 
properties of the MW. Since many of those estimates are 
made using the Johnson-Gousins passband system, we 
will often rely on our results transformed to this system 
in order to make direct comparisons; these are available 
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TABLE 4 

Photometric Properties eor the Milky Way: Rest-frame 2=0.1 .Johnson-Cousins Passbands 


Property 

Corrected Value 

Bias Removed 

9/9M* 



(mag) 

(mag) 

(IQ-IO Mq'^) 

(mag M“° yr) 

- 5 log ft 

-20.10«;“ 

0.252 ± 0.011 

-0.11 

-0.64 

o iMs - 5 log ft 

-19.96t°;« 

0.196 ±0.010 

-0.12 

-0.56 

o iMv - 5 log ft 

-20.47«;" 

0.136 ±0.009 

-0.12 

-0.33 

— 5 log h 

-20.98«;« 

0.143 ±0.008 

-0.12 

-0.25 

o iMj - 5 log ft 

-21.60«;« 

0.139 ±0.008 

-0.13 

-0.36 

O-if^U-V) 

0.604/°; 

0.099 ± 0.0059 

0.02 

-0.10 

B) 

U.Ui^_Q QQO 

0.055 ± 0.0034 

0.01 

-0.01 

- V) 

U.O-^D_q Qg2 

0.037 ±0.0031 

0.01 

-0.09 


U.OiQ_o.o43 

0.010 ±0.0010 

0.00 

-0.04 


U.OO/ _Q Q4y 

0.006 ± 0.0009 

0.01 

0.01 


Note. — Values in this table are determined from analyzing the distributions of properties for Milky Way analogs, but after transforming 
SDSS ugriz measurements to 17BViiL-equivalent values on an object—by-object basis using the kcorrect software. As a reminder, UBVRI 
magnitudes are on the Vega system, whereas ugriz magnitudes are on the AB system. 


TABLE 5 

Global Stellar Mass-to-light Ratios for the Milky Way 


Rest- 

frame 

u 

Y* 

9 

y* 
i r 

'Y'* 

'Y'* 

2 = 0 

i.yu_Q gQ 

J--yo_o g4 

l-OO_0 49 

1 4Q+0-48 

1 11 +0.32 
J-.ll_o.32 

2 = 0.1 

-| yy+1.61 

‘ ‘ -0.83 

1 QQ+O-Sl 

J--y'J_0.68 

to lO 

d d 

+ 1 

GO 

-*-•^^-0.44 

l-^O_0.37 

Rest- 

frame 

'Y'* 

^ u 

'Y'* 

^ B 

T/ 


TJ 

2 = 0 

1.0b_Q gQ 

J--oy_0.65 

^■°^-0.58 

-‘^•'^-‘--0.48 
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Note. — See the end of ij^for details on the calculation of these 
values, which are expressed m units of Mq/Lq. 


in Tables |3|]^ First and foremost, our transformed esti¬ 
mate of°{B — V) = 0.744^0^054 is in excellent agreement 
with the widely used vdK86 measurement of 0.83 ± 0.15, 
consistent at the ~0.5 ct level. Our result indicates a 
slightly bluer color for the MW with a smaller uncer¬ 
tainty by a factor of ~3. 

As mentioned in ^ earlier measurements yielded much 
bluer color estimates for the MW than vdK86, and hence 
also much bluer than the estimate we have presented 
here. The dV&P two-component model produced a color 
estimate oi B — V = 0.53 ± 0.05, which is inconsistent 
at nearly the 3cr level with our estimate. The B&S two- 
component model yields B — V = 0.45; given the lack 
of any error estimates, this is difficult to compare to our 
value, though again significantly bluer. Bahcall (19861 
advised using a ±0.2 mag margin of error when compar- 
ing colors to the model, given the wide variety of system¬ 
atic uncertainties existing in the data at that time. If we 
use this as the error estimate for the B&S model value, we 
find that our result is redder by 0.29 ± 0.21 magnitudes, 
making these estimates inconsistent at the ^lAa signif¬ 
icance level. It is possible that the tension between the 
dV&P and B&S color estimates and the one presented 


here would be reduced if the two-component models em¬ 
ployed were updated to more current constraints on the 
Galaxy’s stellar populations. 

In dV83, an estimated color oi B — V = 0.53 ± 0.04 
is quoted for the MW, obtained by averaging the ob¬ 
served colors of nearby Sb/c type s. This value appears 


to originate from data in Table 4 of de Vaucouleurs (1977 
hereafter dV77), which indicates tnat the distribution o: 
corrected colors for a sample of 70 Sbc galaxies is de¬ 
scribed hy B — V = 0.564±0.066. The measurements for 
each object are tabulated in the Second Reference Cat¬ 


alogu e of Bright Galaxies (RC2; de Vaucouleurs et al. 


19761) , which collected extragalactic data published since 


the l930s. Each B — V color measurement in RC2 is 
corrected to the asymptotic total light from each galaxy 
using a Laplace-Gauss integral technique, as a function 
of its morphological type (T) and the effective aperture 
diameter {Ae) containing 50% of its total light (in some 
cases, the B — V color is transformed from measurements 
in different passbands). Each “total” B — V color is then 
corrected to zero Galactic extinction via a model of the 
Galactic dust as a function of coordinates (^,6), to zero 
internal extinction based on a model of inclination red¬ 
dening as a function of T and isophotal axis ratio, and 
to the z = 0 rest frame via a iF-correction modeled as a 
function of T and z. 

Given the difficulties of these corrections, as well as 
the challenges of properly intercalibrating photographic 
and photoelectric measurements from a wide variety of 
sources, it is likely that there could be significant system¬ 
atic errors in this mean B — V estimate. Eurthermore, 
dV83 assumes T = 4 for the MW (no less than 2.5 and no 
more than 5.5) and quotes the rate of change of the mean 
corrected color along the T sequence near T = 4 to be 
-0.10. While T = 4 (or equivalently Sb/c) fits well with 
the Galactic bulge-to-total ratio of 0.15 we have found in 
LN15, the uncertainty in the MW’s morphological type 
will still represent an additional source of uncertainty 
that appears not to have been included in the error esti¬ 
mate from dV83. We can therefore only treat the uncer¬ 
tainties quoted in the dV83 measurement as a lower limit. 
If we instead consider the value oi B — V = 0.564 ± 0.066 
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from dV77, this is bluer than our nominal result by 0.18 
mag and inconsistent at the ~2(T significance level. We 
believe the tension between our color measurement (or 
that of vdK86) and the estimates from dV77 and dV83 
would be relieved if the sources of uncertainty described 
above were included in this estimate. 


For comparison, Fukugita et al. (19951 performed a 
similar analysis for galaxy types across the Hubble se¬ 
quence by comparing synthetic colors measured from 
galaxy SEDs to broadband photometry taken from the 
Third Reference Cat alogue of Bright Galaxies (de Vau- 
couleurs et al. [1991). Listed in their Table 2, the au¬ 
thors found that the average B — V color for 676 Sb/c 
types (using only objects with |5| > 30°, but applying 
no reddening correction) is 0.68 ± 0.14, which is in ex- 
cellent agreement with our MW result. More recently, 
Fernandez Lorenzo et al. (|2012|) investigated the colors 
ot isolated galaxies m the AMlOA sample that are also 
found in SDSS-DR8. This sample included 466 galaxies, 
two-thirds of which were classihed as Sb/c. Similarly to 
the methods we employ, the authors used model mag¬ 
nitudes that were corrected for Galactic dust extinction 
and AT-corrected to z=0 rest-frame passbands. Listed in 
their Table 3, they found that the median ^(g — r) color 
for Sb/c types is 0.65 ± 0.09, which compares well with 
our MW ^(g — r) estimate of 0.682tQ g 5 g. That table also 
provides colors for a variety of Sb/c galaxy samples; these 
subsets vary in local environment and redshift range, but 
all yield color estimates that agree with our MW value 
at or below the ~ltT significance threshold. 

6.2. Comparisons to Earlier Absolute Magnitude 
Measurements 

Comparisons of absolute magnitudes require more 
care, as they require additional assumptions that are 
prone to systematic error, particularly the value oi h = 
F7o/(100km s“^Mpc“^) used to bring extragalactic dis¬ 
tance estimates (determined from z) and measurements 
based on absolute distances (in pc) onto a common scale. 
For the following discussion we adopt h = 0.7. The 
vdK86 study yielded estimates of Mb = —20.3 ±0.2 and 
(when combined with his B — V estimate) My = — 21.1± 
0.3; these compare well with our slightly brighter esti¬ 
mates of = -20.84^°;^^ and ^My = -21.5lt!^;i^^, 
which are consistent at the ~lcr level. The B&S two- 
component model yields Mb = —20.1 and My = —20.5, 
measurably dimmer than the results we have found, 
though again hard to compare to with no error esti¬ 
mates given. The dV&P two-component model, on the 
other hand, produced Mb = —20.2 ± 0.15 (dV83) and 
B -V = 0.53 ± 0.05, leading to My = -20.7 ± 0.16; 
these values are inconsistent with our ^Mb and re¬ 
sults at the ^1.5 g and ~2cr levels, re spectively. 

More recently, Flynn et al. (2006) analyzed Hippar- 
cos and Tycho data for the local disk and extrapolated 
using an exponential disk model (in combination with 
earlier bulge luminosity estimates) to determine Mj = 
—22.3±0.17. This compares well with our brighter value 
of -22.61", and is consistent with it at the ^0.8cr 


level. Also, Liu et al. (2011) converted the best -to-date 
Vega-calibra ted My measurement for the MW ( [van den 

OOOn I TiTf/~» nn A 1 1 Vvtti o l-vor\l n 0 • J- rr> VkO -n <-1 


Bergh 2000) into an AB-calibrated absolute —band 
magnitude of -21.97 (with no error estimate given); this 


is within ~lcr of our estimate of = —21.55/jg’3g, 

but brighter, rather than fainter. 

6.3. Comparisons to Luminosity Function 
Measurements 


Blanton et al. (2003b) determined luminosity functions 
for galaxies in all SUSS passbands using the SDSS Early 
Data Release. These were determined as the Schechter 
function that fits best to the distribution of Petrosian 
absolute magnitudes of galaxies, converted to the AB 
system, AT-corrected to z = 0.1 rest-frame ugriz pass- 
bands, and corrected for galaxy luminosity evolution; 
this should compare well with the cmodel absolute mag¬ 
nitudes used in this work after correcting them for the lu¬ 
minosity evolution since z = 0.1. The Schechter function 
is parameterized by the characteristic absolute magni¬ 
tude, M* (not to be confused with the total stellar mass 
which we have denoted as M*), which provides a mea¬ 
sure of where the luminosity function transitions from 
being well fit by a power law into an exponential drop¬ 
off. Thus, galaxies with increasing absolute magnitude 
beyond Af* rapidly become more rare. To compare our 
results, we add O.IQ to the results listed in Table of 
this paper, where Q is the appropriate correction in units 
of magnit ude per unit redshift f or each band as listed in 
Table 3 of Blanton et al. (2003b), and then subtract from 
this quantity the appropriate M* value for each band as 
listed in their Table 2. Based on this work, we find that 
the the MW is brighter than their M* by 0.50 ± 0.64, 
0.48 ± 0.48, 0.18 ± 0.37, 0.18 ± 0.38, and 0.15 ± 0.39 
magnitudes in the ^-^ugriz bands, respectively, essen¬ 
tially showing the MW to be consistent with iiA in all 
bands at the 5.1 o' sig nificance level. Similarly, Montero- 


Dorta fc Prad a (2009) reproduced the analysis of Blanton 


et al. 1^03bp using SPSS Data Release 6, which provides 


larger redshift-complete samples of galaxies and incorpo¬ 
rates improved reductions of SDSS imaging data. How¬ 
ever, the luminosity functions that result from this work 
neglect any correction for the evolution of galaxies, as 
its impact is estimated to be very small for the redshift 
ranges used (i.e., z < 0.2). To compare our results with 
this work, we subtract M* listed for the appropriate band 
in their Table 2 from our values in Table of this paper. 
Based on this work, we find that the MW is brighter than 
M, by 1.13 ± 0.63, 0.54 ± 0.48, 0.07 ± 0.37, 0.23 ± 0.38, 
and 0.01 ± 0.39 magnitudes in the ^-^ugriz bands, re¬ 
spectively. Again, we hnd the MW to be consistent with 
their M* in nearly all bands at the ^Icr confidence level, 
and hence is comparable in luminosity to A* galaxies in 
the nearby universe. 


6.4. Comparisons to Green-valley Definitions 

In Figure we have presented the MW’s location in 
the ^{g — r) vs. ^Mr plane, demonstrating that it falls in 
the intermediate region between the blue-cloud and red- 
sequence populations. Our results are consistent with 
the hypothesis posed by Mil that the MW could be a 
member of the “green-valley” population. In Figure 
we have produced an updated °(u — r) vs. logM* dia¬ 
gram modeled on of Figure 1 of Mil, showing the vast 
improvement in our constraints compared to those from 
prior measurements. Here, the MW lies bluer than their 
definition of the green valley, i.e., the region within 0.1 
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log(M*/Mo)-21og(V0.7) 

Fig. 9. — Updated version of Figure 1 from Mil showing the 
Milky Way’s corrected position (red point and purple Icr ellipse) 
in — r) vs. M* space, where again our new constraints are a 
dramatic improvement upon and consistent with the prior measure¬ 
ments (gray dashed—dotted lines) utilized by Mil (updated her e 
to the color transform^ons for galaxies from |Cook et al.||2014|| . 
Comparing with Figure]^ the green valley becomes much stretched 
out in ^(u — r) color space. The dark shaded green region follows 
the same pre scription as Mil, using the empirically derived [Baldry | 
|et al.| (2006 I color division line with a diO.l *^(u — r) offset bec- 
ond, the ligAt shaded green region is the green valley as defined by 
S14 for SDSS galaxies after the effects of dust are removed; this 
provides a suitable comparison for our dust—corrected Milky Way 
results, whereas the grayscale contours do not reflect this correc¬ 
tion. A similar story emerges as from Figure]^ our Galaxy likely 
resides in the saddle of the bimodal color distribution of galaxies 
in the local universe. Measured externally, it would appear redder 
than the majority of spiral galaxies, yet bluer than most ellipticals. 
This makes the Milky Way one of most massive, brightest, and red¬ 
dest of spiral galaxies with appreciable star formation today. 


mag of the SDSS color division line determined bylBaldryl 

et^(|5^. 

However, more recent work by S14, using X-corrected 
and dust-corrected DR7 magnitudes for galaxies at 
0.02 < z < 0.05 (comparable to our sample at 0.03 < z < 
0.09), defines the green valley to be —0.75 < °(u — r) — 
O.251og(M*/M0) < —0.24; this definition would indicate 
that t he MW is i n fact a green-valley galaxy in this dia¬ 
gram. I Jin et al. (2014) present a definition of the green 
valley that avoids dust-reddening effects by using face- 
on nearby galaxies with DR7 magnitudes X-corrected 
to z = 0.1. They define the center of the valley to be 
r) = —0.121(°-^Mr — 51og/i) —0.061 (with no range 
given). Given the uncertainties in our measurements, we 
find that the MW is bluer than this line by 0.25 ± 0.21 
mag, consistent with it at the ~1.2 (t confidence level. 


Mendez et al. (2011) define the green-valley region of 
the — B) vs. Mb plane for AEGIS galaxies to be 
within a ±0.1 mag vertical offset of the line ^{U — B) = 
—0.0189(°Mb — 5 log h) — 0.32 (where we have converted 
from AB to Vega magnitudes). We find that our results 
place the MW redder than this line by 0.089±0.070 mag; 


hence the MW mi ght be considered a gr een-valley galaxy 
by this definition. Willmer et al. (2006) present a similar 
GMD division line tor red and blue galaxies measured 
in the DEEP2 Redshift Survey. They define this line as 
0([/ -B) = -0.032(OMb - 51og(/r/0.7)) - 0.585, where 
we have included small corrections to reflect the AB-to- 
Vega magnitude conversions from kcorrect that have 
been employed in this study. We find that our results 
place the MW redder than this line by 0.067 ± 0.071 
mag. 

It is interesting to note that in the color-magnitude 
plane shown in Figure]^ none of the MW analogs appear 
in the peak of the blue-cloud region where prototypical 
blue, star-forming spirals reside, which would hint that 
our Galaxy, too, very likely does not fit that mold. This 
is contrary, however, to what one finds in the SFR-M* 
plane shown in Figure]^ the vast majority of the MWAS 
lie in the blue cloud or just below. Based on its color, 
if seen from outside, the MW would likely be defined as 
a member of the green valley. In ^[g — r), in fact, it is 
likely very close to the minimum-density region of color 
space. However, based on its M* and M*, it appears to 
fall just off the blue cloud, if it is not actually a member 
of it. It thus provides a cautionary example: objects may 
fall in the green-valley region of parameter space for a 
variety of reasons, especially when only optical (and not 
UV) color is considered. 

6.5. Comparisons to Earlier Mass-to-light Ratio 
Measurements 

In Table we have presented new estimates of the 
global stellar mass-to-light ratio, T*, of the MW in 
SDSS ugriz passbands in the z=0 and 0.1 rest frames, 
as well as these results transformed to Johnson-Gousins 
UBVRI passbands . The most releva nt study we can 
compare these to is Flynn et al. (2006), which presented 
direct estimates of I'"' tor the local Galactic disk by 
accounting for the mass and luminosity budget in the 
“solar cylinder” (i.e., the column of stellar material at 
Rq). This work primarily relied on fitting their Tuorla 
Galactic model to data taken in the Hipparcos and Ty¬ 
cho surveys (reaching out to ^200 pc), which was shown 
to match well with the Heidelberg model-independent 
analysis of the much more shallow (<25-50 pc) Gata- 
logue of Nearby Stars. They found Ty = 1.5 ± 0.2 
M 0 /L 0 , and then used color conversion derived from 
Hipparcos/Tycho data to obtain = 1.4±0.2 M 0 /L 0 
and Tj = 1.2 ± 0.2 M 0 /L 0 . We note that if we were 
to update these to reflect the solar absolute magnitudes 
and colors we have employed herein (~0.03 mag differ¬ 
ences), they would increase by ~3%, well below the Icr 
uncertainties. Regardless, we find that these values com¬ 
pare well with our global MW results of = 1.89tQg5 
M 0 /L 0 , = 1.86;°j9 M 0 /L 0 , and T* = 1.29;°:f7 

M 0 / Lq , whic h are larger than but consistent with the 
corresponding |Flynn et al. estimates at the ^ 0.7cr, 0.6cr, 
and 0.2cr levels, respectively. 

One should keep in mind that, whereas the |Flynn et al.| 
estimates describe the disk itself, our results represent 
the global (disk±bulge) values and hence are expected 
to be larger to some extent, especially in the B- and V- 
bands, as they include the contribution from older stars 
in the Galactic nucleus. We can illustrate this further, 
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and hence make a more apples-to-apples comparison, by 
making the following back-of-the-envelope calculation. 
First, for the subset of our volume-limited sample that 
has /deV > 0.95, which constitutes ^33,000 highly bulge- 
dominated or elliptical galaxies, we find a distribution 
of values that is well approximated as a Gaussian 
described by 4.1 ± 0.9 Mq/Lq (after multiplying by a 
fac tor of 1.5 to convert from Kroupa to Salpeter IMF; 
cf. Fukugita et al. 1998). Second, in LN15 we have de- 
termined the buige-to-total ratio of stellar mass in the 
MW to be B/T = 0.15 ± 0.02. By combi ning our 
estimate for spheroidal components with the [Flynn et al.| 
estimate for the Galactic disk, using the LN15 esti- 
mate of B/T to calculate a mass-weighted average for 
both components, we find a global mass-to-light ratio of 
Ty = 1.81 ±0.19 M-q/Lq. This is in excellent agreement 
with our result and is consistent with it at the ^O.lcr 
level. Doing the analog ous calculatio ns in the V- and 
/-bands, the remaining Flynn et al.j disk values corre¬ 
spond to global values of Ty = 1.80 ± 0.20 Mq/Lq and 
Tj = 1.39 ± 0.18 M 0 /L 0 , which are again in excellent 
agreement with our results, consistent with them at the 
O.lcr and 0.2cr levels, respectively. 


6.6. Conclusions and Future Studies 

Overall, since the vas t majority of spiral galaxies pop¬ 
ulate the b lue cloud (Strateva et al. 2001| Blanton 
eral]|2003a| S14), our results imply that the Mw ranks 


amongst the most luminous, yet reddest of spirals in the 
local universe. Based on a variety of empirical definitions 
in the literature, our results show that it is likely that 
the MW would be classified as a green-valley galaxy if 
viewed from the outside, generally taken to indicate that 
it would be in a transitional evolutionary stage. Again, 
this is contrary to what we find in the SFR-M* plane 
shown in Figure ^ where we hnd that the MW lies very 
near, if not on, the main sequence of star-forming galax¬ 
ies. Apparently, even when the impact of dust effects is 
accounted for, the green valley can be misleading when 
using it to generally characterize the galaxies it contains 
(cf. S14). It is safe to say that our Galaxy’s SFR is in 
a state of decline; the MW produces only ~1.65 solar 
masses of new stars per year, even though it is amongst 
the brightest and most massive of late types. Our find¬ 
ings support the emerging consensus view of the MW; 
one in which it is not the prototypical, blue spiral it was 
commonly thought to be just a decade ago, but is instead 
si milar to the pas sive, red spiral population investigated 


Gortese (2012). In fact, based on the demographics 


of late types presented by S14, if our Galaxy truly lies 
in the green valley then its photometric properties would 
be representative of only ~19% of the spiral galaxies in 
the nearby universe. It is beyond the scope of this paper 
to discuss what evolutionary histories may produce an 
optically red, yet still star-forming MW (or equivalently 
late types tha t appear in the green valley), but we refe r 


the r eader to Hammer et al. ( 2007 ), Yin et al. (2009), 
Mil, Mendez et al. l 201l| , Jin et al. ( 2014[ ), ana S14 tor 


insightful discussions. 

In following papers we will incorporate other well- 


studied MW parameters into our technique of studying 
the MW via its analogs, including morphological type, 
bulge-to-total ratio, and disk scale length. In doing 
so, we should be able to constrain other photometric 
properties that cannot be directly measured (e.g., the 
central surface brightness and global Sersic index), but 
are commonly measured for other galaxies. Additionally, 
we will integrate UV-wavelength data from GALEX to 
more accurately assess whether the MW belongs in the 
green valley (cf. Wyder et al. 2007), as well as utilize 
WISE data to investigate its near-lR properties. The 
same multi-wavelength estimates for our nearest MW- 
like neighbor, M31, the Andromeda Galaxy, whose prox¬ 
imity and thus brightness can cause saturation effects in 
survey data, would also likely improve from the analog 
analysis method. Lastly, the sample of MW analogs we 
have obtained for this work could also be used to explore 
other properties with new observations; e.g., an ancil¬ 
lary target program is now underway to provide integral 
field observations for a subsample of the MWAS studied 


herein as a part of the SDSS-IV MaNGA survey (Bundy 


et al. 2015). In a soon-to-follow paper, we will use the 


results of this paper and those from LN15 to explore the 
position of the MW in a variety of scaling relations for 
disk galaxies. 
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